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Electro-osmotic  drag  and  back  diffusion  are  the  primary  water  transport  mechanisms  in  PEMFC  (Proton 
Exchange  Membrane  Fuel  Cell)  electrolytes.  These  two  phenomena  occur  competitively  in  the  membrane, 
and  ultimately  determine  the  net  water  movement.  The  chemical  compositions  of  reactants  and  product 
(i.e.,  H2,  02,  and  H20)  in  a  porous  catalyst  layer  vary  with  respect  to  the  electro-chemical  reactions  and 
water  transport  through  the  membrane.  The  tendency  of  the  chemical  compositions  was  estimated  by 
analyzing  the  net  water  transport  coefficient  (or),  defined  as  the  ratio  of  the  reaction  rate  to  the  water 
transport  rate.  New  criteria  were  suggested  for  predicting  species  mole  fractions  from  the  flow  direction, 
and  these  were  validated  by  CFD  analysis.  The  hydrogen  mole  fraction  had  different  tendencies  to  rise  or 
fall  based  on  the  flow  direction  at  a  =  0.5,  while  the  oxygen  mole  fraction  extreme  was  located  at  a  =  -0.75. 
Finally,  a  was  shown  to  influence  the  membrane  conductivity  and  activation  losses,  which  are  the  main 
factors  that  contribute  to  fuel  cell  performance. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Worldwide  energy  consumption  is  sharply  rising,  and  new  envi¬ 
ronmental  regulations  will  soon  be  implemented  in  an  effort  to 
reduce  atmospheric  carbon  dioxide  levels,  which  can  influence 
global  warming.  Many  countries  are  concentrating  on  developing 
clean  and  renewable  energy  systems.  Fuel  cells  present  an  environ¬ 
mentally  friendly  and  efficient  energy  conversion  system  that  has 
been  considered  as  an  alternative  power  source.  Additionally,  fuel 
cells  have  a  high  potential  for  commercialization  due  to  continuing 
research.  PEMFCs  (Proton  Exchange  Membrane  Fuel  Cell)  are  low- 
temperature  fuel  cells  that  can  be  used  for  FCVs  (Fuel  Cell  Vehicles) 
or  RPG  (Residential  Power  Generation)  systems.  It  is  proposed  that 
FCVs  will  be  the  biggest  market  in  10  years;  thus,  major  automobile 
companies  throughout  the  world  have  concentrated  on  commer¬ 
cializing  PEMFCs.  Flowever,  long-term  stability,  sub-zero  operation, 
and  competitive  prices  must  be  achieved. 

To  maintain  stable  and  efficient  operation  of  PEMFCs,  irre¬ 
versible  losses  such  as  activation  losses,  ohmic  losses,  and 
concentration  losses,  must  be  minimized.  Activation  losses  are 
unavoidable  when  generating  current  from  a  fuel  cell,  and  these  are 
affected  by  temperature,  pressure,  catalyst  materials,  and  the  mole 
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fractions  of  reactants.  The  mole  fractions  of  reactants  are  related 
to  the  mass  transport  in  a  bipolar  plate  (BPP)  channel  and  a  gas 
diffusion  layer  (GDL).  Since  concentration  losses  originate  from 
the  process  of  supplying  reactants  and  removing  products,  con¬ 
centration  losses  are  partially  related  to  activation  losses  [1].  It  is 
known  that  concentration  losses  account  for  a  significant  portion 
of  the  overall  irreversible  losses,  and  thus,  a  uniform  distribution  of 
reactant  must  be  achieved  to  improve  the  performance  of  PEMFCs 
[2-7]. 

Mass  transport  in  the  cathode  is  more  complex  than  in  the 
anode,  due  to  the  complexity  of  oxygen  diffusion  and  water 
removal.  Additionally,  cathode  activation  losses  are  generally  larger 
than  those  in  the  anode  because  the  barrier  energy  for  the  oxygen 
reduction  reaction  (ORR)  is  higher  than  for  the  hydrogen  oxidation 
reaction  (FIOR).  Therefore,  a  uniform  oxygen  distribution  is  impor¬ 
tant  for  reducing  activation  losses  and  concentration  losses,  and  has 
been  investigated  extensively  [8-16]. 

Humidity  is  also  an  important  parameter  in  reducing  ohmic 
losses  because  the  resistance  of  the  membrane  is  determined  by 
the  relative  humidity.  Springer  et  al.  measured  and  modeled  the  ion 
conductivity  of  a  117  Nafion  membrane  with  varying  humidity  [17]. 
Their  results  have  been  widely  employed  for  estimating  membrane 
conductivities.  Buchi  and  Scherer  compared  the  humidity  effect 
at  both  electrodes.  Two  and  four  membranes  were  laid  to  detect 
the  resistance  of  each  membrane.  Finally,  the  authors  found  that 
humidity  at  the  anode  is  more  influential  than  at  the  cathode  when 


D.  Lee,  J.  Bae  /  Journal  of  Power  Sources  191  (2009)  390-399 


391 


Nomenclature 

ak  water  activity  in  stream  k 

ACv  specific  surface  area  of  control  volume  (nrr 1 ) 

Areacv  surface  area  of  control  volume  (m2 ) 

Cn  concentration  of  species  n  (mol  nrr3 ) 

Dnj  diffusion  coefficient  of  species  n  in  mixture  j 
(m2  s-1) 

D'n  j  diffusion  coefficient  of  species  n  in  mixture  j  through 

porous  media  (m2  s-1 ) 
d  diffusion  adjustment  (%) 

DH2  j  diffusion  coefficient  of  H2  in  liquid-water  film 
(6.3e-9m2  s-1) 

Dq2  \  diffusion  coefficient  of  02  in  liquid-water  film 
(2.4e-9m2  s-1) 

F  Faraday  constant  (96,487  C  mol-1 ) 

Hh2,i  Henry’s  law  constant  for  H2  in  liquid- water  film 
(8.9e9Pa) 

Hq2  j  Henry’s  law  constant  for  02  in  liquid-water  film 
(2.12elOPa) 

hrxn  enthalpy  for  water  formation  (kj  kmol-1 ) 
hfg  enthalpy  for  water  vaporization  (kj  kg-1 ) 

/  local  current  density  (A  nr 2 ) 

/0,n  exchange  current  density  for  n  (A  m-2 ) 
mnk  mass  fraction  of  species  n  in  stream  k 
Mn  equivalent  weight  of  n  (kg  mol-1 ) 
massn  mass  of  species  n  (kg) 
rid  electro-osmotic  drag  coefficient 

P  pressure  (Pa) 

Pn  partial  pressure  of  species  n  (Pa) 

R  universal  gas  constant  (8.314J  mol-1 I<-1 ) 

r  condensation  rate  (nr3  s-1 ) 

S  source  terms 

She  heat  source  term  from  electrochemical  losses 

Shp  heat  source  term  from  water  phase  change 

£m  membrane  thickness  (m) 

f/./<  liquid-water  film  on  the  k  (anode  or  cathode)  (m) 

T  temperature  (I<) 

u  velocity  (ms-1) 

V0c  open-circuit  voltage  (V) 

Vceii  cell  voltage  (V) 

Xnk  mole  fraction  of  species  n  in  stream  k 

[n]  molar  flux  of  species  n  (mol  s-1  m-2) 

k  thermal  conductivity  (W  nr 1  K-1 ) 

Greek  symbols 
£  porosity 

r\  kinetic  overpotential  (V) 

/a  dynamic  viscosity  (kg  s-1  nr 1 ) 

p  density  (kg nr3) 

a  net  water  flux  per  proton  flux  through  membrane 

ak  kinetic  transfer  coefficient  for  reaction  k 
ex m  membrane  conductivity  (S  nr 1 ) 

permeability  (m) 

Superscript 

sat  saturation 

Subscript 
a  anode 

c  cathode 

e  electrochemical  reaction 

p  phase  change 

H2  hydrogen 


02  oxygen 

1  liquid 

v  water  vapor 

H20  water 

§  dummy  variable  for  direction  x,  y,  or  z 

rxn  reaction 

oxi  oxidation  at  anode  by  reaction 
red  reduction  at  cathode  by  reaction 

drag  electro-osmotic  drag  through  membrane  from 
anode  to  cathode 

back  back  diffusion  through  membrane  from  cathode  to 
anode 

gen  water  generation  at  the  cathode  by  reaction 

van  vanished  water  at  cathode 


the  membrane  thickness  is  less  than  50p,m  [18].  In  other  words, 
hydrogen  has  to  be  well  humidified  to  minimize  ohmic  losses. 

As  discussed,  a  uniform  oxygen  distribution  at  the  cathode  and 
high  relative  humidity  at  the  anode  are  necessary  conditions  to 
improve  the  performance  of  PEMFCs.  These  two  conditions  can  be 
satisfied  without  external  humidification  when  the  membrane  is 
self-humidified  by  water  generated  from  fuel-cell  reactions.  Back 
diffusion  can  decrease  activation  losses  and  increase  membrane 
conductivity.  Vengatesan  et  al.  suggested  that  the  performance  of 
PEMFCs  can  be  enhanced  by  using  a  thin  cast  membrane  [19].  In 
this  setup,  the  membrane  conductivity  can  be  properly  maintained 
by  the  back  diffusion.  The  possibility  of  internal  humidification  was 
also  mentioned  by  Buchi  and  Srinivasan  [20].  It  was  also  reported 
that  back  diffusion  is  intensified  by  adding  hydrophilic  Si02  parti¬ 
cles  to  the  catalyst  layers  at  anode,  and  cell  performance  is  improved 
at  low  relative  humidity  [21].  However,  the  relationship  between 
water  movement  through  the  membrane  and  cell  performance  has 
not  been  clearly  defined  in  the  literatures. 

The  net  water  movement  through  the  membrane  and  the  fuel¬ 
cell  reactions  are  linked  to  each  other  and  can  cause  the  chemical 
composition  to  vary.  For  example,  the  composition  of  the  anode 
stream  is  determined  by  the  water  transport  from  electro-osmotic 
drag  and  the  HOR.  In  contrast,  species  mole  fractions  in  the  cath¬ 
ode  can  be  estimated  from  the  water  back-diffusion  rate  and  the 
water-generation  rate  (or  alternatively,  the  oxygen  reaction  rate). 
Consequently,  in  this  study,  each  species  mole  fraction  in  the  anode 
and/or  the  cathode  and  the  net  water  movement  are  shown  to  be 
related  and  have  the  influence  for  the  membrane  conductivity  and 
activation  losses.  Criteria  are  suggested  for  predicting  the  rise  and 
fall  in  the  oxygen  mole  fraction,  the  water  mole  fraction  in  the 
anode,  and  the  irreversible  losses.  CFD  analysis  was  used  to  validate 
the  relationship  between  the  mole  fraction  change  of  each  species 
in  the  catalyst  layer  and  the  net  water  flux  through  the  membrane. 


2.  Mathematical  formulation 

Mass,  momentum,  and  energy  conservation  equations  were 
used  as  the  governing  equations  for  simulating  a  PEMFC.  Species 
transport  equations  for  hydrogen,  oxygen,  and  water  were  also  used 
to  predict  the  changes  in  mole  fractions  in  the  anode  and  cath¬ 
ode.  The  following  assumptions  were  applied  in  the  PEMFC  model: 
steady  state,  laminar  flow  in  the  channel,  high  electronic  conduc¬ 
tivity  in  the  bipolar  plate  (BPP)  for  uniform  galvanic  potential,  and 
homogeneous  two-phase  flow  [8]. 

The  mass  conservation  equations  are  as  follows: 


V-(p2)  =  Sm 


(1) 
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5m  —  5h2  +  SWVp  +  5wip  +  Sawve  «3tZ  —  Z3 

5m  —  5o2  +  5wvp  +  5wip  +  5C  wve  atz  =  Z2 

where, 

(1.1) 

(1.2) 

5H2  -  2p  ^h2A:v  atz-z3 

(2.1) 

So2  =  -  I-^Lp-Mo2A„  at z  =  z2 

(2.2) 

a(x  v) 

Sawve  =  p  I{x,y)MH2oAcv  atz  =  z3 

(2.3) 

_  1+2  a(x,y)„  ... 

5cwve  —  2  p  h2oA;v  atz  —  z2 

(2.4) 

5wip  —  — 5WVp 

MH20^nofy(masSnofv/Mnofv)  ( P^_  —  P wv  \ 

(hm  V  ^  / 


The  mass  change  of  each  species  due  to  reaction  is  described  by 
Eqs.  (2.1  )-(2.4).  The  phase  change  between  water  vapor  and  liquid 
water  is  affected  by  the  saturation  pressure  and  defined  by  Eq.  (2.5), 
where  r  is  the  condensation  rate  (1  m-3  s-1 ).  The  liquid  water  can 
move  to  the  vapor  phase  by  evaporation  if  the  partial  pressure  of 
liquid  water  is  higher  than  the  saturation  pressure. 

The  momentum  conservation  equations  are  as  follows: 


V  •  ( puu )  =  -  VP  +  V  •  (/zVu)  +  Sp, 

(3) 

where, 

Sp  =  - ^ i/  atz2  <  z  <  zv  andzi  <  z  <  z3 

(3.1) 

The  momentum  sink  terms,  defined  by  Eq.  (3.1),  originate  from 
Darcy’s  law  and  are  used  for  modeling  the  pressure  drop  in 
the  porous  catalyst  layer.  The  permeability  {ft)  is  very  small 
(1.0e-12m2)  and  assumed  to  be  isotropic.  Therefore,  the  source 
terms  are  derived  from  the  Navier-Stokes  equation  because  the 
other  terms  become  negligible  [5]. 

The  energy  conservation  equations  are  as  follows: 

V (puh)  =  V  •  (kVT)  +  Shp  +  5he  (4) 

where, 

c  U  fl{x,y)xAcv\ 

She  =  hrxn  X  ( - — - 

-  (J(x, y)  X  Vceii  X  ACV)Z=Z2  at z  =  0  (4.1 ) 


Shp  —  *^wlp  x  ^fg*  (4*2) 

The  difference  between  the  total  chemical  energy  of  the  reac¬ 
tants  and  the  electrical  work  was  assumed  to  be  released  as  heat. 
Since  /irxn  is  for  the  overall  chemical  reaction  between  H2  and  02, 
the  heat  of  reaction  is  assumed  to  be  generated  at  the  middle  of  the 
membrane,  and  not  on  the  catalyst  surface.  Latent  heat  due  to  phase 
change  is  also  contained  in  the  heat  source  terms  and  is  evaluated 
by  the  enthalpy  of  formation  of  water,  hfg.  The  flow  and  thermal 
properties  in  the  GDL,  the  BPP,  and  the  membrane  are  shown  in 
Table  1. 

The  species  transport  equations  are  as  follows: 

V  •  (pmnu)  =  V  •  C/n)  +  Sn,  (5) 

where  n  denotes  the  species,  H2,  02,  water  vapor,  or  liquid  water. 
The  change  in  the  mass  fraction  resulting  in  electrochemical  reac¬ 
tions  is  reflected  by  the  source  terms  in  the  species  equations.  These 
source  terms  are  the  same  as  for  the  mass  conservation  equations. 


Table  1 

Flow  and  thermal  parameter  in  the  GDL,  the  BPP,  and  the  membrane  [23]. 


Porosity  in  the  GDL 
Permeability  in  the  GDL 
Diffusion  adjustment 
GDL  thermal  conductivity 
BPP  thermal  conductivity 
Membrane  thermal  conductivity 


0.7 

1.0e-12  m 
50% 

0.21  Wm-1  K"1 
5.7  Wm-1  K"1 
0.15  Wm-1  K"1 


The  equations  and  parameters  used  to  model  the  electrochem¬ 
ical  phenomena  are  shown  in  Tables  2  and  3.  The  water  activity 
indicates  the  humidity  at  the  membrane  (Eq.  (8)),  and  the  water 
content  (A.)  represents  the  number  of  water  molecules  per  sulfuric 
acid  bonds  in  the  membrane  (Eq.  (9)).  The  membrane  conductivity 
was  determined  by  the  water  activity  and  water  content  (Eq.  (14)). 
It  is  known  that  ohmic  losses  are  more  significantly  affected  by  the 
humidity  at  the  anode  than  at  the  cathode  when  the  membrane  is 
less  than  50  jxm  [18].  Therefore,  the  membrane  conductivity  (crm) 
can  be  evaluated  from  the  water  content,  based  on  the  anode  water 
activity. 

Two  main  water  transport  phenomena  occur  in  the  electrolyte 
of  PEMFCs.  Electro-osmotic  drag  (n^)  results  from  proton  move¬ 
ments  that  promote  water  transport  from  the  anode  to  the  cathode 
(Eq.  (10)).  Back  diffusion  is  caused  by  a  concentration  difference 
that  promotes  water  diffusion  from  the  cathode  to  the  anode,  and 
was  modeled  by  the  water  diffusivity  of  the  membrane  (Eq.  (11)) 
and  the  water  concentration  (Eq.  (12)).  The  net  water  flux  per  pro¬ 
ton  through  the  membrane  (a)  is  defined  by  the  flux  difference  of 
electro-osmotic  drag  and  back  diffusion  (Eq.  (13)).  By  definition, 
water  is  transported  from  the  anode  to  the  cathode  if  a  is  positive, 
and  vice  versa.  All  of  the  source  terms  in  Eq.  (2.1)-(2.5)  are  linked 
because  a  is  included  in  each. 

The  activation  loss  ( r /)  was  estimated  by  the  Tafel  equation,  as 
shown  in  Eq.  (15)  [1  ].  The  exchange  current  density  (J0)  is  modified 
by  the  hydrogen  or  oxygen  mole  fraction  in  the  natural  logarithm  of 
Tafel  equation.  The  concentration  loss  can  be  included  in  the  acti¬ 
vation  loss  by  modifying  the  exchange  current  density.  The  cathode 
activation  loss  {r]c)is  generally  higher  than  in  the  anode  due  to  the 
complexity  of  the  ORR.  Therefore,  the  total  activation  loss  is  deter¬ 
mined  dominantly  by  the  cathode  activation  loss  and  the  oxygen 
mole  fraction. 

Current  densities  were  calculated  from  Eq.  (16),  considering  the 
open-circuit  voltage  (OCV)  and  irreversible  losses.  The  local  mem¬ 
brane  conductivity  and  activation  losses  are  the  main  influencing 
factors  for  the  current  densities. 

Liquid  water  is  formed  when  the  water  partial  pressure  is  higher 
than  the  saturation  pressure.  It  is  assumed  that  the  catalyst  layer  is 
partially  covered  by  a  thin  liquid  film  [23].  Hydrogen  and  oxygen 
dissolve  in  the  film  according  to  Henry’s  law  (Eq.  (17)).  The  diffusion 
length  of  soluble  gas  is  determined  by  the  thickness  of  the  liquid- 
water  film.  The  thickness  of  liquid  water  is  modeled  by  Eq.  (18),  as 
discussed  in  the  appendix  of  Lee  et  al.  [16]. 

All  governing  equations  and  appropriate  boundary  conditions 
were  solved  using  the  commercial  CFD  tool  STAR  CD  ver.  3.26,  with 
custom  user  subroutines. 


3.  New  criteria  for  prediction  of  the  species  mole  fractions 

The  species  mole  fractions  can  be  estimated  in  a  porous 
catalyst  layer  by  considering  the  reaction  rate  and  the  water  trans¬ 
port  through  the  membrane.  That  is,  the  variation  of  species 
mole  fractions  can  be  predicted  by  the  net  water  flux  per  pro¬ 
ton  (a)  because  the  number  of  protons  related  to  the  reaction 
rate  and  net  water  movements  through  membrane  are  contained 
in  oi. 


Table  2 

Equations  for  the  electrochemical  reactions. 
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Variables 


Equations 


Diffusion  mass  flux 

Binary  diffusion  coefficient  [22,23] 

Water  activity 

Water  content  in  the  membrane  [17] 
Electro-osmotic  drag  coefficient  [11  ] 

Water  diffusion  coefficient  in  the  membrane  [17] 


Water  concentration  on  the  surfaces  of  the  MEA  [1  ] 
Net  water  flux  per  proton  [23] 

Local  membrane  conductivity  [13] 

Local  overpotential  [11] 

Local  current  density  [11  ] 

Gas  solubility  in  the  liquid-water  film  [8] 

Water  film  thickness  [8] 


h,n  =  -PDt,n 


dmK,n 

% 


PDnJ(x,y ) 


=  3.64  x  10“ 


(Pc,n  X  PcJ)1/3(Tc,n  X  TcJ)5/12(l/M„  +  1/Mj)1/2 
D'n  j(x,  y)  =  d  x  Dnj(x,  y )  in  porous  media 

Xw  ,KP(x,  y) 

aK=—j^t - 

w  ,K 

X  —  0.043  +  17.81aa  -  39.85a;;  +  36a3  (0  <  aa  <  1) 
X  =  14  +  1.4(afl  -  1)  (1  <  aa  <  3) 
nd  =  0.0029A2  +  0.05A  -  3.4  x  10“19 


Ux,y) 

\/Tc,nTc, 


Dw=Dxx  exp  (  2416  (  — 
10 


{2m(rn-rh)) 


A  <2 


Dx=  ■ 


10_1°(1  +  2(A  -  2)),  2  <  A  <  3 

10_1°(3  -  1 .67(A  -  3)),  3  <  X  <  4.5 


1.25  x  10“ 


4.5  <  A 


/(X?y)  =  am[X,y)(V0C  -  VcELL  -  Tj(x,  y)) 


J(x,y) 

2F 

K*,y) 


Mh2  —  Ph2 
Mq2  =  / 


f/,k  = 


H2’’  l - - ) 

/  Xo2  ,zglifi c  PzgU, c  /Ho2 , 1  —  Xo2,z2  \ 
4F  —  2=Po2Do2,i  ^ - - - j 

(^mass  „) 


(6) 

34 

(7) 

(B) 

(9) 

(10) 

(11) 


Cw,Jf(x,y)=  xAaorc 

Mm,  dry 

<*(*,  y)  =  na(x,  y)  -  w(x,  y)  C"'4^r)  -  Cwia(x.y) 

<Tm(x,  y)  =  (0.514A  -  0.326)  x  exp  ( 1268  (  ^  -  )  ) 

n(x,  y)  =  r\c(x,  y)  +  iu(x,y)  =  In  f-  /(^y)  )  +  In  f  J(^y)  )  (15) 

0Ccr  \w,02*02,z2  J  \iO,H2AH2,z3  / 


(12) 

(13) 

(14) 


epwiAreaa 


(16) 

(17) 

(18) 


Six  variables  were  used  to  compare  the  conversion  rate  of  each 
species  (Fig.  1 ).  Hydrogen  is  oxidized  and  two  protons  are  generated 
at  the  anode  catalyst  layer.  Oxygen  is  reduced  with  protons  which 
are  dragged  from  the  anode,  and  converted  into  water  molecules  at 
the  cathode  catalyst  layer.  [H2]0Xi  is  the  oxidation  rate  of  hydrogen, 
and  [02  ]red  is  the  reduction  rates  of  oxygen.  Water  is  dragged  from 
the  anode  to  the  cathode  ([H20]drag),  or  in  the  opposite  direction 
([H20]back  =  —  [H20]drag).  [H20]gen  is  the  water-generation  rate  at 
the  cathode,  and  [H+]drag  is  the  electro-osmotic  drag  rate  of  protons 
through  the  membrane  from  the  anode  to  the  cathode. 

Kinetic  equations  were  derived  (Eqs.  (19.1)  and  (20.1))  from  the 
electrochemical  reactions  at  the  anode  (19)  and  the  cathode  (20).  a 
was  also  derived  (Eqs.  (19.2)  and  (20.2))  at  both  catalyst  layers  from 


Table  3 

Parameters  for  the  electrochemical  reaction  [23]. 


Dry  membrane  density  (pm,dry) 

2  kg  ur 3 

Equivalent  weight  of  dry  membrane  (Mmdry) 

1.1  kg  mol-1 

Exchange  current  density  for  hydrogen  (/0,h2  ) 

0.2  A  cm-2 

Exchange  current  density  for  oxygen  (/0,o2 ) 

0.02  A  cm-2 

Anode  transfer  coefficient  (aa) 

1.2 

Cathode  transfer  coefficient  (ac) 

0.6 

the  drag  rate  of  protons  and  the  water  transport  rates. 


H2  2H+  +  2e  at  anode  catalyst  layer 

(19) 

1H  1  IH+ldrag 

1*^2  Joxi  —  2 

(19.1) 

[H20]drag 
“  [H+]drag 

(19.2) 

2H+  +  2e“  +  jC>2  4*  H2O  at  cathode  catalyst  layer 

(20) 

[H+]drag  In  ,  [H20]gen 

4  =  1^2]  red  =  2 

(20.1) 

[H20]back 

[H+]drag 

(20.2) 

New  criteria  for  predicting  the  species  mole  fractions  are  sug¬ 
gested  by  Eqs.  (21)  and  (22).  These  formulations  are  derived  from 
the  combination  of  kinetic  equations  and  the  definition  of  a.  These 
criteria  will  be  referred  to  as  the  ‘alpha  criteria’  throughout  this 
paper.  If  a  >  0.5,  [H20]drag  is  greater  than  [H2]oxi  at  the  anode  cat¬ 
alyst  layer  in  Eq.  (21),  and  vice  versa.  If  -0.5  <  a  <0,  [H20]gen  is 
greater  than  [H20]back  at  the  cathode  catalyst  layer  in  Eq.  (22). 


[H2U  j  [HjO]** 

1  Anode  1 

[Hl^g  ; 

- 1 —  | 

|H*W 

|  Cathode  | 

[o2U 

[H20]back 


[FhOjgcn 


Fig.  1.  Conceptual  diagram  of  the  electrochemical  reaction  within  the  catalyst  layer  and  water  movement  through  the  membrane. 
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Fig.  2.  Schematic  diagram  of  the  alpha  criteria:  five  regions  of  a  (a),  mole  fraction 
changes  at  the  anode  (b),  and  mole  fraction  changes  at  the  cathode  (c). 

[H20]Van,  defined  by  the  difference  between  [H20]gen  and  [H20]back, 
is  smaller  than  [02]red  when  a  <-0.75.  The  direction  of  water 
movement  through  the  membrane  is  determined  at  o'  =  0  by  def¬ 
inition.  Therefore,  a  can  be  separated  into  five  distinct  regions  by 
these  four  points  (o' =  0.5,  0,  -0.5,  -0.75)  to  provide  comparative 
rates  at  each  electrode. 

[H20]drag  =  (2o')[H2]oxi  (21) 

[H20]van  =  [H20]back  -  [H20]gen  =  —2(2a  4-  l)[02]red  (22) 

A  schematic  diagram  of  the  alpha  criteria  is  shown  in  Fig.  2. 
The  five  distinct  regions,  regions  (i)-(v),  are  shown  in  Fig.  2(a).  The 
change  in  species  mole  fractions  in  the  anode  stream  is  represented 
in  Fig.  2(b).  The  tendency  of  species  composition  in  the  anode 
stream  can  be  characterized  at  a  =  0.5  since  [H20]Van  is  comparable 
t°  [H2]0Xj  at  o' =  0.5  from  Eq.  (21).  In  the  region  where  o' >0.5,  the 
water  mole  fraction  may  decrease,  or  the  hydrogen  mole  fraction 
may  increase  through  the  anode  stream  because  [H2  ]oxi  <  [H20]drag. 
If  o'  <  0.5,  the  anode  water  mole  fraction  may  increase,  or  the  hydro¬ 
gen  mole  fraction  may  decrease.  The  rise  and  fall  of  the  membrane 
conductivity  will  be  influenced  by  these  tendencies  of  the  water 
mole  fraction  in  anode.  The  change  in  species  mole  fractions  in 
the  cathode  stream  is  represented  in  Fig.  2(c).  The  tendency  of 
species  composition  in  the  cathode  stream  can  be  characterized  at 
o' =  -0.75  since  [H20]drag  is  comparable  to  [02]red  at  o' =  -0.75  from 
Eq.  (22).  The  water  mole  fraction  may  increase  or  the  oxygen  mole 
fraction  may  decrease  through  the  cathode  stream  in  the  region 
where  a  >  -0.75  because  [FI20]Van  <  [02  ]red.  If  a  <  -0.75,  the  cath¬ 
ode  water  mole  fraction  may  decrease,  or  the  oxygen  mole  fraction 
may  increase.  The  rise  and  fall  of  activation  losses  will  be  influ¬ 
enced  by  these  tendencies  of  oxygen  mole  fraction  in  the  cathode. 
In  summary,  changes  in  the  behavior  of  each  species  along  the  flow 
direction  are  dissimilar  before  and  after  a  =  0.5  for  the  anode,  and 
before  and  after  a  =  -0.75  for  the  cathode.  These  tendencies  can 
be  linked  to  membrane  conductivity,  anode  and  cathode  activation 
losses,  and  current  density.  Among  these  parameters,  membrane 


Table  4 

Boundary  conditions. 

Inlet  flow  temperature 

70  °C 

Inlet  humidity  (case  number) 

Anode:  100% 

Cathode:  70%  (1),  80%  (2),  90%  (3) 

Inlet  humidity  (case  number) 

Anode:  70%  (4),  75%  (5),  80%  (6) 
Cathode:  100% 

Stoichiometry 

6  (anode  and  cathode) 

Temperature  at  BPP  surface 

70  °C 

conductivity  and  the  cathode  activation  losses  will  be  investigated 
using  the  alpha  criteria. 

4.  Computational  domain 

4.1.  Grids  without  gas  diffusion  layers  (GDLs) 

The  computational  domain  for  validating  the  alpha  criteria  is 
shown  in  Fig.  3.  The  mid-plane  of  the  membrane  is  located  at  z  =  0. 
The  membrane  thickness  is  50  p,m,  which  is  the  distance  between 
Z\  and  Zv.  The  thicknesses  of  the  catalyst  layer  and  the  chan¬ 
nel  are  100  and  250  pirn,  respectively.  The  interfaces  between  the 
channel  and  the  porous  catalyst  layer  are  at  Z2  and  Z3.  The  flow 
channel  is  approximately  2  cm  long  in  the  x  direction.  A  total  of 
20,800  nodes  were  used  to  model  the  PEMFC  (120  x  8  x  5  nodes 
for  both  the  anode  and  the  cathode  channels  and  100  x  8  x  4  nodes 
(membrane),  3  nodes  (both  BPPs),  or  2  nodes  (both  catalyst  lay- 
ers)). 

To  accurately  predict  the  behavior  of  the  species  mole  fractions 
using  our  new  criteria,  gas  diffusion  layers  are  not  introduced  in 
Fig.  3.  When  the  gas  diffusion  layers  are  included,  concentration 
gradients  appear  at  the  porous  medium.  It  was  found  that  this 
results  in  an  error  in  the  estimation  of  the  alpha  criteria.  Addi¬ 
tionally,  the  ribs  of  the  interconnector  are  not  included  in  the 
computational  domain  to  consider  the  one-dimensional  changes 
of  species  mole  fractions  in  the  x  direction.  The  stoichiometry  at 
both  the  anode  and  the  cathode  channels  was  set  at  6.  High  sto¬ 
ichiometries  and  the  small  volume  of  the  porous  medium  (from 
neglecting  the  GDLs)  made  the  species  concentration  uniform  in 
the  z  direction.  These  limitations  will  be  addressed  in  Section  5 
with  an  analysis  of  the  results  including  the  GDLs  and  the  ribs  under 
realistic  stoichiometric  conditions. 

Hydrogen  and  water,  provided  by  a  humidifier,  were  injected  in 
the  anode  channel.  Oxygen,  instead  of  air,  and  steam  were  injected 
in  the  cathode  channel.  The  humidity  conditions  were  controlled 
to  determine  the  net  water  movement  through  the  membrane.  Six 
cases  were  applied  to  the  grids  without  GDLs  to  verify  the  alpha 
criteria  (Table  4).  Electro-osmotic  drag  will  be  more  active  than  back 
diffusion  in  cases  1  -3  due  to  the  high  relative  humidity  at  the  anode. 
However,  water  will  move  from  the  cathode  to  the  anode  by  back 
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Fig.  3.  Computational  domain  of  a  single  straight  channel  flow  field  for  a  PEMFC  without  GDLs:  (a)  geometry  and  (b)  cross-sectional  view. 
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Fig.  4.  Computational  domain  of  a  single  straight  channel  flow  field  for  a  PEMFC  with  GDLs:  (a)  geometry  and  (b)  cross-sectional  view. 


diffusion  in  cases  4-6  because  of  the  high  relative  humidity  at  the 
cathode.  The  flow  directions  in  the  anode  and  cathode  were  the 
same  in  the  positive  x  direction.  The  current  density  was  0.3  A  cm-2 
for  all  cases.  Results  were  extracted  from  the  lines  located  aty  =  0 
andz=Z!  orZr. 


4.2.  Grids  for  a  real  PEMFC  with  GDLs 

Negligible  vertical  flow  between  the  porous  catalyst  layer 
and  the  channel  and  one-dimensional  concentration  change  were 
assumed,  in  order  to  limit  the  computational  domain.  However,  in 
real  PEMFCs,  the  channel  thickness  is  usually  greater  than  1  mm, 
and  the  GDL  is  located  between  the  channel  and  the  catalyst  layer. 
Therefore,  reactant  and  product  concentrations  are  non-uniform 
throughout  the  diffusion  medium.  Additionally,  species  concentra¬ 
tion  changes  are  no  longer  one-dimensional  if  the  ribs  are  included 
in  the  computational  domain. 

A  new  PEMFC  grid  with  GDLs  and  ribs  was  made  for  comparison 
with  Fig.  3.  The  grid  and  dimensions,  which  has  32,352  nodes,  are 
shown  in  Fig.  4.  The  channel  length  was  4  cm.  The  humidity  condi¬ 
tions  were  the  same  as  case  1  of  Table  4.  More  realistic  stoichiome¬ 
tries  were  used,  1.2  in  the  anode  and  2.0  in  the  cathode.  Humidified 
hydrogen  and  air  were  injected  in  the  anode  and  cathode  channel. 
The  current  density  was  0.6  A  cm-2  and  the  results  were  extracted 
from  the  center  line  in  the  x  direction  at  the  catalyst  layer. 

5.  Results  and  discussion 


mole  fraction  at  the  anode  increased,  except  in  region  (i),  as  previ¬ 
ously  predicted  by  the  alpha  criteria.  The  increasing  tendency  of  the 
anode  water  mole  fraction  along  the  flow  direction  was  maintained 
in  the  outlet  region  in  Fig.  6(a)  and  over  the  entire  area  in  Fig.  6(b). 
This  is  because  the  HOR  rate  is  higher  than  the  water  transport  rate, 
or  water  is  accumulated  by  back  diffusion. 

The  membrane  conductivity  showed  the  same  tendency  as  the 
water  mole  fraction  at  the  anode.  The  conductivity  in  the  inlet 
region  decreased  due  to  the  drop  in  the  anode  water  mole  frac¬ 
tion,  in  Fig.  7(a).  However,  the  membrane  conductivity  increased 


(a) 


x/L 


5.1.  Net  water  flux  per  proton  (a) 

The  net  water  transfer  coefficient  per  proton  is  shown  in  Fig.  5. 
Water  moved  from  the  anode  to  the  cathode  at  the  inlet  region,  as 
shown  in  Fig.  5(a),  due  to  the  fully  saturated  humidity  condition  at 
anode.  Water  diffused  from  the  cathode  to  the  anode  at  the  outlet 
region.  Back  diffusion  was  very  active  over  the  entire  area  in  cases 
4-6  (Fig.  5(b)).  Dashed  lines  are  drawn  at  a  =  0.5  and  -0.75,  where 
the  water  mole  fraction  at  the  anode  and  the  oxygen  mole  fraction 
at  the  cathode  show  different  tendencies  because  of  the  alpha  crite¬ 
ria.  The  locations  of  a  =  0.5  were  at  x/L  =  0.28  for  case  1, 0.23  for  case 
2,  and  0.155  for  case  3,  and  the  points  in  which  a  =  -0.75  were  rep¬ 
resented  at  x/L  =  0.29  for  case  4, 0.235  for  case  5,  and  0.14  for  case  6. 

5.2.  Water  mole  fraction  at  the  anode  and  the  membrane 
conductivity 

The  water  mole  fraction  at  the  anode  is  shown  in  Fig.  6.  There  is 
a  decreasing  area  of  water  mole  fraction  in  Fig.  6(a)  because  water 
is  removed  relatively  faster  than  hydrogen  in  region  (i).  The  water 


<b>  x/L 


Fig.  5.  Net  water  flux  per  proton  along  the  x  axis  at  the  catalyst  layers:  cases  1-3  (a) 
and  cases  4-6  (b). 
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Fig.  6.  Water  mole  fraction  along  the  x  axis  at  the  anode  catalyst  layers:  cases  1-3 
(a)  and  cases  4-6  (b). 


in  the  outlet  region  in  Fig.  7(a)  and  over  the  entire  area  in 
Fig.  7(b). 

5.3.  Oxygen  mole  fraction  and  cathode  activation  losses 

The  oxygen  mole  fraction  at  the  cathode  is  shown  in  Fig.  8.  The 
oxygen  mole  fraction  decreased  in  Fig.  8(a)  and  (b),  except  in  region 
(v).  In  these  cases,  the  rate  of  water  accumulation  from  electro- 
osmotic  drag  or  back  diffusion  was  lower  than  the  ORR  rate.  In 
contrast,  the  oxygen  mole  fraction  increased  in  region  (v),  where 
strong  back  diffusion  occurred  (Fig.  8(b)). 

The  tendency  of  the  oxygen  mole  fraction  has  an  influence  on 
the  cathode  activation  losses.  Cathode  activation  losses  increased, 
as  shown  in  Fig.  9(a)  and  (b),  except  in  region  (v).  This  was  due  to 
the  diminution  of  the  oxygen  mole  fraction.  However,  the  cathode 
activation  loss  increased  in  region  (v)  (Fig.  9(b))  as  opposed  to  being 
reduced  due  to  the  rise  in  the  oxygen  mole  fraction. 

From  Eq.  (15),  it  is  observed  that  the  current  density  and  the 
oxygen  mole  fraction  are  the  main  factors  influencing  the  activation 
losses.  The  activation  loss  will  increase  due  to  the  increase  in  current 
density  or  the  decrease  in  the  oxygen  mole  fraction.  In  cases  1-3, 
the  current  density  weakened  along  the  flow  direction  (Fig.  10(a)). 
The  tendencies  of  current  density  were  expected  to  cause  the  cath¬ 
ode  activation  losses  to  decrease;  however,  the  cathode  activation 
losses  increased  in  Fig.  9(a).  Therefore,  it  can  be  concluded  that  the 
oxygen  mole  fraction  played  a  dominant  role  in  cases  1-3.  In  con¬ 
trast,  the  current  density  increased  along  the  flow  direction  in  cases 
4-6  (Fig.  10(b)),  and  cathode  activation  losses  increased  in  region 
(v),  in  Fig.  8(b).  Therefore,  current  density  proved  to  be  more  influ¬ 


ential  in  cathode  activation  losses  than  the  oxygen  mole  fraction 
in  the  inlet  region,  shown  in  Fig.  8(b).  It  can  be  deduced  that  the 
cathode  activation  losses  must  be  evaluated  with  respect  to  both 
oxygen  mole  fraction  and  current  density. 

5.4.  General  rules  for  the  ‘alpha  criteria’ 

It  was  observed  that  the  species  mole  fractions  in  the  catalyst 
layer  are  influenced  by  the  net  water  transfer  coefficient  per  proton. 
General  rules  for  the  ‘alpha  criteria’  can  be  summarized  as  follow: 

Region  (i)  (0.5  <  a): 

Anode.  0<  [H2]oxi  ^  [F^OJdrag 

Cathode:  [H20]van  =  [H2(3]back  —  [H20]gen  <  0  <  [C^Jred 

The  water  transport  rate  from  the  anode  to  the  cathode  is  higher 
than  the  HOR  rate  at  the  anode  catalyst  layer.  It  seems  that  water 
is  removed  relatively  faster  than  hydrogen  in  the  anode  stream  of 
the  BPP  channel  or  the  GDLs.  Therefore,  the  hydrogen  mole  frac¬ 
tion  increases  and  the  water  mole  fraction  decreases  at  the  anode 
along  the  flow  direction  in  this  region.  The  membrane  conductiv¬ 
ity  may  weaken  due  to  this  tendency  of  the  anode  water  mole 
fraction. 

In  the  cathode,  water  is  significantly  accumulated  with  genera¬ 
tion  and  movement  from  the  anode.  Therefore,  as  the  water  mole 
fraction  increases,  the  oxygen  mole  fraction  decreases  in  this  region. 
The  decrease  in  the  oxygen  mole  fraction  affects  the  increase  in  the 
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Fig.  8.  Oxygen  mole  fraction  along  the  x  axis  at  the  cathode  catalyst  layers:  cases 
1-3  (a)  and  cases  4-6  (b). 

cathode  activation  losses.  Finally,  the  total  activation  losses  may 
increase  in  this  region. 


Region  (ii)  (0  <  a  <  0.5): 

Anode.  0<  [H20](jrag  <  [ H2 ]0xi 

Cathode:  [FbOJvan  =  [H20]back  -  [H20]gen  <0<  [O2 ]red 

The  reaction  rate  of  hydrogen  is  higher  than  the  water  transport 
rate  from  the  anode  to  the  cathode.  Therefore,  hydrogen  is  removed 
relatively  faster  than  water  in  the  anode  stream.  The  hydrogen 
mole  fraction  decreases  and  the  water  mole  fraction  at  the  anode 
increases  along  the  flow  direction  in  this  region.  The  membrane 
conductivity  may  be  enhanced  due  to  the  increase  in  the  water 
mole  fraction  in  the  anode. 

In  the  cathode,  water  is  also  accumulated  with  generation  and 
movement  from  the  anode  since  water  is  transported  from  the 
anode  to  the  cathode.  Therefore,  the  tendencies  of  the  oxygen  mole 
fraction  and  the  water  mole  fraction  are  similar  to  those  in  region 
(i).  The  oxygen  mole  fraction  decreases  and  the  water  mole  frac¬ 
tion  increases.  The  decrease  in  the  oxygen  mole  fraction  affects  the 
increase  in  the  cathode  activation  losses.  Thus,  the  total  activation 
losses  may  also  increase  in  this  region. 


Region  (iii)-(v)  at  the  anode  [a  <  0): 

Anode.  [H2O]cjrag<0'<  [H2]oxi 

When  a  is  negative,  water  is  transported  from  the  cathode  to 
the  anode.  Therefore,  water  is  accumulated  by  back  diffusion  at  the 


anode.  Hydrogen  is  oxidized  and  removed  at  the  catalyst  surface, 
but  water  is  accumulated  in  the  anode  stream.  Therefore,  the  ten¬ 
dency  of  the  water  and  hydrogen  mole  fractions  for  negative  a  is 
similar  to  those  in  region  (ii),  at  the  anode.  The  hydrogen  mole  frac¬ 
tion  decreases,  and  the  water  mole  fraction  at  the  anode  increases 
along  the  flow  direction.  The  membrane  conductivity  may  increase 
due  to  the  behavior  of  the  water  mole  fraction  at  the  anode  when 
a<0. 

Region  (iii)  at  the  cathode  (-0.5  <  a  <  0): 

Cathode:  [H20]Van  =  [H20]back  -  [H20]gen  <0<  [O2 ]red 

In  region  (iii),  the  water-generation  rate  is  greater  than 
the  back-diffusion  rate  at  the  cathode.  Water  is  accumulated 
at  the  cathode  for  this  reason,  although  water  movement 
through  the  membrane  occurs  from  the  cathode  to  the  anode. 
The  oxygen  mole  fraction  decreases  along  the  flow  direction 
at  the  cathode  because  the  oxygen  is  reduced  at  the  cat¬ 
alyst.  Thus,  the  total  activation  loss  may  increase  in  region 
(iii). 

Region  (iv)  at  the  cathode  (-0.75  <  a  <  -0.5): 

Cathode:  0<  [H20]van  =  [H20]t>ack  -  [H20]gen  <  [C^lred 

In  this  region,  the  back-diffusion  rate  is  greater  than  the  water- 
generation  rate.  That  is,  water  is  removed  at  the  cathode  catalyst 
layer  by  back  diffusion  even  though  water  is  generated  from  the 
ORR.  However,  oxygen  is  more  quickly  reduced  than  water  is 
removed.  For  this  reason,  the  oxygen  mole  fraction  decreases,  and 
the  water  mole  fraction  increases  along  the  flow  direction  at  the 


Fig.  9.  Cathode  activation  losses  along  the  x  axis:  cases  1-3  (a)  and  cases  4-6  (b). 
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cathode.  With  the  change  in  the  oxygen  mole  fraction,  the  total 
activation  losses  may  still  increase  in  region  (iv). 

Region  (v)  at  the  cathode  ( a  <  -0.75): 

Cathode:  0<  [02]red  <  [H2O]  van  “  [H20]baci<  —  [H20]gen 

Back  diffusion  is  so  strong  in  this  region  that  water  is  more 
quickly  removed  at  the  cathode  catalyst  layer  than  oxygen.  The  oxy¬ 
gen  mole  fraction  increases  in  this  region  while  the  water  mole 
fraction  decreases  along  the  flow  direction.  These  tendencies  may 
affect  the  diminution  of  the  total  activation  loss. 

5.5.  The  relationship  between  water  movement  through  the 
membrane  and  species  mole  fraction  in  a  real  PEMFC  with  GDLs 
and  BPP  ribs  at  low  stoichiometric  conditions 

The  results  shown  in  Figs.  5-10  describe  the  situation  for  PEMFCs 
operating  without  GDL  at  high  stoichiometry.  Flowever,  GDLs  are 
generally  sandwiched  between  the  channel  and  the  porous  catalyst 
layer,  and  the  stoichiometry  is  not  greater  than  4  in  a  real  PEMFC. 
Therefore,  the  relationship  between  water  movement  through  the 
membrane  and  the  species  mole  fractions  at  the  catalyst  layer 
should  be  analyzed  at  realistic  PEMFC  operating  conditions.  To 
revaluate  the  alpha  criteria  for  these  conditions,  a  new  grid  was 
used,  as  shown  in  Fig.  4.  As  previously  mentioned  in  Section  4.2, 
the  relative  humidity  at  the  inlet  was  the  same  as  case  1,  and  the 
current  density  was  0.6  A  cm-2.  The  stoichiometries  were  set  at  1.2 
in  the  anode  and  2.0  in  the  cathode. 

The  net  water  transfer  flux  per  proton  (a)  is  shown  in  Fig.  11(a). 
The  point  a  =  0.5  was  located  at  x/L  =  0.07  in  the  PEMFC  with  GDLs 
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Fig.  11.  Comparison  between  case  1  and  the  new  case  with  GDL  and  low  stoichiom¬ 
etry:  net  water  transfer  coefficient  per  proton  (a)  (a)  and  hydrogen  mole  fraction  in 
the  anode  catalyst  layer  (b). 


and  ribs.  The  maximum  hydrogen  mole  fraction  was  at  x/L  =  0.15,  as 
shown  in  Fig.  11(b).  Some  errors  occurred  in  predicting  the  maxi¬ 
mum  hydrogen  mole  fraction  from  the  alpha  criteria.  As  GDLs  were 
added  in  a  new  computational  domain  at  low  stoichiometry,  the 
concentration  gradient  through  the  diffusion  layer  became  influ¬ 
ential  for  the  species  composition  in  the  catalyst  layer.  The  species 
mole  fractions  in  the  catalyst  layer  were  not  determined  only  by 
proton  and  water  movement  through  the  membrane,  but  also  by  the 
diffusion  rate  of  hydrogen  and  oxygen  from  the  GDL.  For  this  reason, 
it  is  difficult  to  accurately  estimate  the  tendency  of  species  mole 
fractions  at  the  catalyst  layer  by  considering  only  alpha  criteria,  as 
shown  in  Fig.  11(b). 

Although  the  alpha  criteria  are  inaccurate  due  to  diffusion 
through  the  GDL  in  a  realistic  PEMFC,  the  comparison  between  the 
water  transport  rate  and  the  reaction  rate  can  be  made  reliably. 
For  the  anode,  the  water  transport  and  reaction  rate  are  derived 
in  Eqs.  (23)  and  (24)  from  the  source  terms  in  the  mass  conser¬ 
vation  equation,  Sh2  and  SaWe.  The  water  transport  rate  and  the 
FIOR  rate  are  shown  in  Fig.  12.  The  molar  flux  of  water  was  greater 
than  that  of  hydrogen  when  a  was  over  0.5,  regardless  of  the  GDLs, 
stoichiometry,  or  geometry,  as  shown  in  Fig.  12(a)  and  (b). 

[H2]oxi  =  -%^  (23) 
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Fig.  10.  Current  density  along  the  x  axis:  cases  1-3  (a)  and  cases  4-6  (b). 
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Fig.  12.  Hydrogen  and  water  molar  flux  through  the  membrane:  case  1  (a)  and  the 
new  grid  with  GDLs  (b). 

6.  Conclusion 

The  influence  of  proton  and  water  transport  through  PEMFC 
membranes  on  the  species  mole  fractions,  activation  losses,  and 
membrane  conductivity  was  investigated.  Criteria  for  forecasting 
the  species  mole  fractions  in  the  catalyst  layer  were  validated  using 
CFD  analysis  when  the  GDL  and  BPP  ribs  were  excluded  from  the 
computational  domain  at  high  stoichiometry.  The  hydrogen  and 
water  mole  fractions  at  the  anode  had  their  extrema  at  a  =  0.5.  The 
water  transport  rate  was  larger  than  the  hydrogen  oxidation  rate 
when  a  was  greater  than  0.5.  Therefore,  the  membrane  conductiv¬ 
ity  increase  when  a  was  below  0.5,  and  vice  versa.  The  oxygen  and 
water  mole  fractions  at  the  cathode  had  their  extrema  at  a  =  -0.75. 
The  oxygen  reduction  rate  was  greater  than  water  for  a  >-0.75. 
Although  the  oxygen  mole  fraction  had  a  significant  effect  on  the 
activation  losses,  the  current  density  also  proved  to  be  a  deter¬ 
mining  parameter.  Both  the  oxygen  mole  fraction  and  the  current 
density  must  be  considered  to  evaluate  the  activation  losses. 

The  prediction  of  species  mole  fractions  from  alpha  criteria  was 
not  exact  in  the  computational  domain  with  GDLs  at  low  stoichiom¬ 
etry.  When  the  GDL  was  sandwiched  between  the  catalyst  layer  and 


the  bipolar  plate,  the  concentration  gradient  through  the  diffusion 
layer  cannot  be  neglected.  That  is,  the  species  mole  fractions  in  the 
catalyst  layer  were  not  determined  only  by  proton  and  water  move¬ 
ment  through  membrane  but  also  by  the  diffusion  rate  of  hydrogen 
and  oxygen  from  the  GDL.  Although  the  alpha  criteria  had  some 
errors  predicting  the  maximum  and  minimum  species  mole  frac¬ 
tions,  comparison  of  the  reduction  rate  of  reactants  and  water  on 
the  membrane  and  the  catalyst  surfaces  was  accurately  estimated 
from  the  alpha  criteria. 

The  diffusion  effect  in  GDLs  was  not  fully  investigated  in  this 
paper.  However,  the  species  mole  fractions  can  be  predicted  more 
accurately  in  the  real  PEMFC  channel  by  alpha  criteria  if  the  diffu¬ 
sion  effect  in  GDLs  is  complemented.  Further  study  will  be  needed 
to  define  constraints,  by  applying  the  alpha  criteria  for  actual  con¬ 
ditions  and  understand  the  relationship  between  the  species  mole 
fractions  and  diffusion  in  the  GDL. 
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